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ABSTRACT 

We present the results of a high resolution optical spectroscopic study of nine white dwarf 
candidate members of Praesepe undertaken with the VLT and UVES. We find, contrary to a 
number of previous studies, that WD0836+201 (LB390, EG59) and WD0837+199 (LB393, 
EG61) are magnetic and non-magnetic white dwarfs respectively. Subsequently, we determine 
the radial velocities for the eight non-magnetic degenerates and provide compelling evidence 
that WD0837+185 is a radial velocity variable and possibly a double-degenerate system. We 
also find that our result for WD0837+218, in conjunction with its projected spatial location 
and position in initial mass-final mass space, argues it is more likely to be a field star than a 
cluster member. After eliminating these two white dwarfs, and WD0836+199 which has no 
clean SDSS photometry, we use the remaining 5 stars to substantiate modern theoretical mass- 
radius relations for white dwarfs. In light of our new results we re-examine the white dwarf 
members of Praesepe and use them to further constrain the initial mass-final mass relation. 
We find a a near monotonic IFMR, which can still be adequately represented by simple linear 
function with only one outlier which may have formed from a blue straggler star. 
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1 INTRODUCTION 

The initial mass-final mass relation (IFMR) is a theoretically pre- 
dicted positive correlation between the main sequence mass of a 
star with M^IO M© and the ma ss of the white dwarf rem nant left 
behind after it has expired (e.g. Ilben Ir. & Renzinill 19831) . Under- 
standing the form of this relation is important since it provides a 
handle on the total amount of gas enriched with He, N and other 
metals that low or intermediate mass stars, which account for 95 
per cent of all stars, return to the interstellar medium at the end of 
their lifecycles. Moreover, the form of the upper end of the IFMR 
is relevant to studies of Type II supernovae as it can provide a con- 
straint on the minimum mass of star that will experience this fate. 

The form of the IFMR is extremely difficult to predict from 
theory alone due to the many complex processes occurring during 
the final phases of stellar evolution (e.g. th ird dredge-up, thermal 
pulses, mass loss; Ilben Jr. & Ren zini 1983). This means that ro- 
bust empirical data are essential for constraining its form. How- 
ever, these are by no means simple to obtain, a significant dif- 
ficulty being the determination of the main sequence mass of 



* Based on observations made with ESO telescopes at the La Silla Paranal 
Observatory under programme ID 076.D-0751 
f E-mail: slc25@star.le.ac.uk 



a star that has long since ceased to exist. This difficulty can 
be a lleviated by using white dwarf member s of open star clus- 
ters teergeron et aljll995t IClaver et alj|200ll; Ipobbie et alj|2004 



Williams et al.l 12004; iKalirai et alj |2005t IDobbie et al.l l2006dlbl ; 
KaliraietalJl2007l;IWilliams &Boltell2007l) to define the IFMR 
(e.glWeidemannl 1 1971 l200dl; IFerrario et al.l 120051 ; IDobbie et al] 
2006a; [Williams & Boltd200llKarirai et alj2008h . Here, since the 
age of the population c an be determine d fr om the location of the 
main sequence turn-off dSandagd 1 19561 ; e.g. lKing&Schulerll2005l 
in Ursa Major) the lifetime, and ultimately the mass, of the progen- 
itor star of any degenerate member can be estimated by calculating 
the difference between the cooling time of the white dwarf and the 
cluster age. 

Nonetheless, until relatively recently, rather fe w white dwarf 
membe rs of open star clusters had been identified. iDobbie et al.l 
(2006a) suggest ~30 associated with 14 clusters, while the IFMR 
derived by IFerrario et alj d2005l) used 40 DA white dwarfs from 7 
open star clusters, including NGC2099. The uncertainties in mem- 
bership status, cluster ages and the relatively large distances of 
these scarce, faint objects resulted in large scatter in the IFMR 
delayer et al. I l200ll ; IFerrario et alJ fcoOS). To begin to address this 
issue and gauge the level of intrinsic scatter in the IFMR, we re- 
cently investigated the white dwarf members of the moderately rich 
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nearby Praesepe open cluster. At a distan ce of 177lo°.; 3 pc (as de - 
termined from Hipparcos measurements, iMermilliod et alj[l997l) . 
it is one of the closest star clusters. I t is slightly me tal rich with 
respect to the Sun ([Fe/H] = +0.11, lAn et alj[2007h . Indeed, as 
both the metallicity and the kinematics of Praesepe are similar to 
those of the Hyades, the former is often touted as a member of the 
Hyades moving group and therefore is ass umed have an age com- 
parable to the latter, r=625±50 Myr (e.g. Iciaver et al.ll200lh . We 
note that this age for the Hyades was derived by comparing model 
isochrones generated from slightly metal enhanced (Z=0.024) stel- 
lar models which included moderate convective overshooting to the 
colours and magnitudes of a sa mple of cluster membe rs selected 
using Hipp arcos astrometric data |Perrvman et aljl998t) . 

In Dobbie et al 1 d2004 [20063T we increased the number of 
probable white dwarf members of Praesepe from five to eleven. 
All thes e candidate members share the proper motion of the clus- 
ter (e.g. Iciaver et al.ll200ll ; lDobbie et alj|2004h . We used estimates 
from the literature and our own measurements of the effective tem- 
perature and the surface gravity of these white dwarfs to determine 
their masses and the masses of their progenitor stars so that we 
could study their locations in initial mass-final mass space. Com- 
bining these with data from a number of other open clusters with 
comparatively well constrained ages and the Sirius binary system, 
we concluded that most stars follow rather closely a monotoni- 
cally increasing IFMR, which can be approximately described by 
a simple linear function over the initial mass range 2.7 Mq to 6 
Mq. Nevertheless, there existed a small number of outliers, all at- 
tributable to the Praesepe cluster. 

Here we present the results of a high resolution optical spec- 
troscopic study of nine out of eleven candidate Praesepe degener- 
ates. In the next section we describe the data acquisition and anal- 
ysis, including the measurement of line core velocity shifts and 
the redetermination of effective temperatures and surface gravi- 
ties. We then utilise the five most appropriate Praesepe degener- 
ates to re-examine the theoretical white dwarf mass-radius relation. 
Subsequently we determine the radial velocities for all eight non- 
magnetic white dwarfs and use these to confirm their membership 
status. Finally, we examine our new results in the context of the 
IFMR. 



2 ACQUISITION AND REDUCTION OF DATA 

High resolution spectroscopy of nine of the eleven DA white dwarf 
candidate members of Praesepe was obtained using the Ultravio- 
let and Visual Echelle Spectrograph (UVES) on the UT2 facility 
of the European Southern Observatory's (ESO) Very Large Tele- 
scope at Cerro Paranal, Chile. The data were acquired in service 
mode within the period 2005/11/16 and 2006/01/19 and using 45 
minute exposure times. As only relaxed constraints were specified 
for the sky conditions these data were generally obtained through 
thin cirrus and in mediocre seeing (^1.4"). 

The observations utilised the UVES parallel mode, where 
cross dispersers CD2 and CD3, in conjunction with dichroicl 
(390+580) and the HER5 and SHP700 filters provided simultane- 
ous wavelength coverage over 370-500nm and 420-680nm in the 
blue and red arms respectively. The 1.5" slit resulted in a nominal 
spectral resolution of A/AA~20000. 

The data were reduced using the UVES pipeline 
which is based on the M IDAS software supplied by ESO 
dBallester & Hensbergej [l995l) . The pipeline was run in OPTIMAL 
mode and an additional custom routine written by R. Napiwotzki 



applied to reduce the magnitude of the ripples which are sometimes 
apparent in these spectra. Calibration frames were taken every 
night science data was acquired. These include format-check 
frames, bias frames, flat fields, an order definition template, and a 
ThAr arc lamp spectrum. In addition to these calibrations, a spec- 
trally featureless DC white dwarf (WD0000-345) was observed 
with the same instrument set-up on 2006/07/03. The spectrum 
of this star was used to remove remaining UVES instrumental 
signature from the data. The completed spectra, which each have 
a S/N~25 per resolution element, can be seen in figure [T] These 
spectra have been median filtered, smoothed and co-added for the 
purpose of this figure. 



2.1 Data analysis 

Our first step in analysing these echelle data was to remove the sig- 
nature of telluric water vapour from around Ha using an absorp- 
tion template constructed from numerous UVES spectra of white 
dwarfs o btained at Paranal Obs ervatory as part of the SpY pro- 
gramme dNapiwotzki et alfcOOll) . Next, the data were compared to 
the predictions of white dwarf model atmospheres using the s pec- 
tral fitting programme FITSB2 (v2.04; lNapiwotzki et ai]|2004l) . 

The grid of pure-H model spectra was calculated using 
the plane-parallel, hydrostatic, non-local thermody namic equilib- 
rium (non-LTE) atmo sphere code TLUSTY, v200 faubenvlfl988l : 



iHubenv & Lanzj|l995 ) and the spectral synthesis code SYNSPEC 
v48 dHubenv & LanzfeooTh . TLUSTY assumes plane-parallel ge- 
ometry and hydrostatic equilibrium. The models include a treat- 
ment for co nvective energy trans port according to the ML2 pre- 



scription of Ber geron et alj d!992t) . adopting a mixing length pa- 
rameter, a=0.6. These calculations utilised a model H-atom which 
incorporates explicitly the eight lowest energy levels and represents 
levels n=9 to 80 by a single superlevel. The dissolution of the high 
lying levels was treated by mean s of th e occupation probability 
formalism of iHummer & Mihalasl J 1988h genera lised to the non- 
LTE atmosphere situation bv lHubenv et al .I d 1994). All calculations 
include the bound-free and free-free opacities of the H~ ion and 
incorporate a full treatment for the blanketing effects of HI lines 
and the Ly man —a, —f3 and —7 satellite opacities as computed by 
N. Allard dAllard et alj|2004l) . During the calculation of the model 
structure the lines of the Lyman and B aimer series were treated by 
means of an approximate Stark profile but in the spectral synthesis 
step detailed profiles for t he B aimer lines were calculated from the 
Stark broadening tables of Lemk^ dl997l) . The grid of model spec- 
tra covered the T eS range of 13000-20000 K in steps of 1000 K and 
log g between 7.5 and 8.5 in steps of 0.1 dex. 



2.2 Determination of effective temperature and gravity 

Each white dwarf has been observed at least three times. We used 
FITSB2 to fit our grid of model spectra to the eight Balmer ab- 
sorption lines ranging from Ha to H10 in each exposure. So as to 
match the instrument resolution, the models were convolved with 
a Gaussian with a full width half maximum of 0.2 A. In addition, 
points in the observed data lying more than 3<r from the model 
were clipped from subsequent iterations of the fitting process. Both 
the starting values and the new values derived here for T e £f and 
log g for each white dwarf are given in Table Q] The errors given 
in Table Q] for T a g and log g are formal fitting errors and are un- 
realistically small as they neglect systematic uncertainties e.g. flat 
fielding errors and model shortcomings. In subsequent discussion 
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Figure 1. The full median filtered, smoothed, stacked UVES spectra for the nine white dwarfs included in this study. Note the Zeeman shifted components of 
the Balmer lines in the spectrum of WD0836+201 (see Figure[3]for more detail). The asterisks mark artifacts caused by dead pixels in the array. There are also 
two breaks in the data of width ~80A at 4580 and 5640A caused by the edges of the CCDs. 



Table 1. ID, effective temperature, lo g g, mass and cooling time for the Praesepe w hite dwarfs with UVES spectroscopy. The previous effective temperature 
and surface gravity estimates are from Dobbie et al. 1 2006a) and lClaver et alj d200ll) . where the values shown are as listed in this literature. The FITSB2 data 
are shown with the formal fitting errors. An uncertainty of 2.3 per cent and 0.07 dex are more realistic for T c ff and log g. Additionally, for Praesepe members, 
we list our estimate of the progenitor mass. The magnetic white dwarf WD0836+201 has been omitted as our models do not include a treatment for magnetism. 
Additionally, no initial mass has been given for WD0837+2 1 8 as it is believed more likely to be a field star than a cluster member. 
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here we follow iNapiwotzki et al.l dl999h and assume an uncertainty 
of 2.3 per cent in T c s and 0.07dex in log g. 

2.3 Determination of line core velocity shifts 

Using our new determinations for the values of T e g and log g (with 
the exception of WD0836+201) as input parameters, a model grid 
which was more finely sampled in wavelength space than the first 
(0.05 A) was used in conjunction with FITSB2 to measure the ve- 
locity shift of the line cores of the Ha and H/3 lines in the spectra of 
each white dwarf. Crude estimates for the line core velocity shifts 
were obtained initially, taking into account the heliocentric correc- 
tion. The line fitting was then re-run using these as starting values 
for the final fit, to minimise the errors and ensure that each result 
was robust. We made reliable estimates of the uncertainties on the 



line shift measuremen ts based on t he the bootstrapping method of 
statistical resampling ( lEfronll 19821) . We note that a detailed study 
of the wavelengths of prominent water absorption lines in a large 
number of UVES spectra obtained as part of the SpY project, which 
uses an identical instrumental set-up to this study, indicates that er- 
ror in the external calibration of the wavelength scale of UVES 
from spectrum to spectrum is at the 0.7 kms -1 level. As an inde- 
pendent check of our results we also performed these fits using a 
different set of sy nthetic spectra (the LTE model atmospheres of 
iKoester et al J200lb . Fortunately, no systematic discrepancies were 
found between the results derived using the two sets of models and 
the differences for any object were well within the measurement 
uncertainties. The line core velocity measurements for each obser- 
vation may be found in Table|2] 

In order to describe the statistical likelihood that a given tar- 
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Table 2. ID, line core shift measurements for each data set (RVi,RVs, RV3, RV4, x 2 ,and logio probability of the white dwarf being a radial velocity variable, 
for the eight non-magnetic Praesepe white dwarfs investigated. The errors take into account the 0.7 kms -1 uncertainty in the UVES wavelength calibration. 
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WD0833+194 


77.13±2.09 
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WD0837+199 
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WD0840+200 


73.76±3.27 
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0.2153261 
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WD0843+184 


93.40±3.08 


86.10±1.99 


87.66±3.10 


87.80±2.06 


4.008115 


-0.584 



get's velocity is constant or otherwise between each observation, 
we calculated a \ 2 statistic for each star based on the prior assump- 
tion that the velocities are stable with time i.e. any object with a \ 2 
below 10 can be considered to have a constant radial velocity. It 
can be seen from Table[2]that while the measurements for most ob- 
jects are consistent with a constant velocity, there is evidence that 
the wavelengths of the line cores in the spectra of WD0837+185 
change from exposure to exposure. We note that in the SpY project 
where velocity measurements have been obtained for a sample of 
1000 white dwarfs, objects with a probability of less than 10 -3 , 
corresponding to one false "detection" due to chance statistical 
fluctuations, are flagged as potential variables. Given our assess- 
ment of the probability of obtaining the tabled radial measurements 
for WD0837+185 by chance «10~ 4 ), we are led to conclude that 
this white dwarf is a radial velocity variable. 



3 RESULTS 

3.1 A possible double degenerate? 

A cursory glance at Figures [2] and [3] appears to substantiate the 
above conclusion that WD0837+185 is a radial velocity vari- 
able. However, we find no compelling direct evidence of any 
cool companion within either the UVES or the Sloan Digital 
Sky Survey (SDSS) spectrum of WD0837+185, the latter ex- 
tending to 9200A. There is also no detection of this ob ject in 
the 2MASS Point Source Catalogue dSkrutskie e t al. 2006]), al- 
though it is listed in the UKIRT Infrared S ky Survey Large 
Area Survey (UKIDSS LAS lWarren et alj|2007t) , with magnitudes 
Z=18.29±0.03, Y=18.49±0.04 and #=18.50±0.14 (on a Vega 
system). Taken in conjunction with the SDSS optical photome- 
try (tj=18.35±0.03, #=18.02±0.02, r=18.32±0.01, i=18.56±0.02, 
z=18.76±0.04, on approximately an AB magnitude system), these 
also provide no evidence for the presence of a cool c ompanion (see 
Figure [Hi. Indeed, based on the DUSTY models of IChabrier et al.l 
( 2000), we set a limit of M coo i ^0.05 M Q , at the distance and age of 
this cluster. Given the r arity of brown dwarf + white dwarf binaries 
(e.g. iFarihi et al .1 2005 estimate the fraction of L dwarf compan- 
ions to white dwarfs to be at <0.5 per cent), this suggests that any 
companion is more likely to be another white dwarf. 

However, there is no evidence in the UVES spectrum which 
would indicate the presence of a He-atmosphere DB companion 
and the observed Balmer line profiles are relatively well matched 
by a synthetic spectrum for a single DA white dwarf. Moreover, 
at the spectroscopically determined effective temperature and sur- 
face gravity, the SDSS magnitudes of WD0837+185 are consistent 
with those of a single degenerate at the distance of Praesepe. The 
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Figure 2. The Ha line profile in the three sets of spectra observed for 
WD0837+185. The red lines are the FITSB2 model fits to the data. The 
core of the lines can be seen moving. These observations were taken on 
2005/1 1/16 (UT1) and 2005/12/14 (UT2 + UT3). The observations are plot- 
ted in chronological order with the earliest observation at the top of the plot. 



models o f iHolberg & Bergeronl $2O0(h indicate an isolated DA of 
this effective temperature and surface gravity to have M 9 ~l 1.6. If 
we were to assume that WD0837+185 had a degenerate companion 
with M 9 Ssl2.5, the absolute magnitude of this hypothetical binary 
would be M 9 Jsll.2. Since the observed magnitude of this object 
is <?=18.02±0.02, this would infer a double degenerate lying well 
beyond the back of the cluster (diJ220 pc). So, if this is a binary, 
the likelihood of Praesepe membership, as indicated by the proper 
motion and the mean radial velocity, would appear to favour a pu- 
tative white dwarf companion which is significantly less luminous 
thanWD0837+185. 

For a cooling time ^625 Myr (i.e. the age of the clus- 
ter) the available data suggests that any de generate companion 
must be massive (M^1.2 M Q ; M„<q2.5: Ulthaus et all 120071: 
IHolberg & Bergeron| [2006), with a low intrinsic luminosity due to 
its small radius. If this is the case, the total mass of this system 
would be larger than the Chandrasekhar limit (~ 1.4 Mq) and this 
binary has the potential to be a progeni t or for a Type la super- 
nova dlben&Tutukovlll984l ; |Livioll2000l: iNapiwotzki et alJl200lL 
120031) . Further study of this system is warranted to better constrain 
the masses of the components, and to determine if the system is 
close enough that it will merge within a Hubble time. We note that 
WD0837+185 sits on the IFMR defined by the bulk of the other 
objects (see Figure [7), so there is nothing in terms of the mass of 
this white dwarf to suggest that this system has experienced close 
binary evolution. 
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Figure 3. As for Figureff] but for the H/3 line profile from the three sets of 
spectra observed for WD0837+185. 
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Figure 5. The full median filtered, smoothed, stacked UVES spectra for 
WD0837+199 and WD0836+201 covering the H/3 line profile. The Zeeman 
shifted components can be seen for WD0836+201. 
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Figure 4. Sloan u, g, r and i, and UKIDSS Z, Y and H band photometry 
of WD0837 + 185, compared to a synthetic DA white dwarf spectrum for 
T e ff = 15076K and log g = 8.306. There is no evidence for excess emis- 
sion due to a cool companion. At an age of 625Myr these data suggest that 
such an object, if present, must have a mass < 0.05Mq . 



3.2 Magnetic white dwarf 

We note here that the identities of WD0837+199 and WD0836+201 
have been confused in som e p revious studies, e.g. IClaver et al.l 
d200lh : IDobbie et alj < |2004T) and lDobbie et al.1 fc006ah . The UVES 
spectrum of WD0836+201 (LB 390, EG59) clearly shows Ze eman 
shifted components (see Figure [5}. This is consistent with iReia 
dl996h who notes his fit to the Ha line core in the spectrum of 
WD0836+201 (LB390) is particularly poor, as would be expected 
if trying to fit Zeeman shifte d com ponents with non-magnetic syn- 
thetic spectra. IClaver et alj d200ll) a lso finds a magn etic white 
dwarf in Praesepe, and in Figure 8 of dClaver et alJl2001D . presents 
the fit to the spectrum, but labels it as EG61 (i.e. WD0837+199, 
LB393) for which th ey determine T eff =1 7,098 K and log #=8.32. 
Moreover, Table 3 of IClaver et al.1 ( 120010 erroneously label EG61 
as WD0836+201. We believe from our data that the m agnetic white 
dwarf is EG59 and th i s objec t has been mislabelled in|Claver et al.1 
d200lh . iDobbie et al.1 d2006j) and iDobbie et al.1 d2004D never~ot> 
served WD0836+201 or WD0837+ 199 but refer to EG 61 as the 
magnetic white dwarf as identified bv lClaver et al. The Zee- 

man shifted Balmer lines in the UVES spectrum of WD0836+201 
due to the substantial magnetic field (B~3 MG), prevents us from 
obtaining a meaningful result via fitting the absorption lines with 
our non-magnetic models. 



3.3 Constraints on the white dwarf mass-radius relation 

The white dwarf mass-radius relation is founded on a Nobel 
Prize winning theoretical de scription of the equat ion of state of 
an electron degenerate gas dChandrasekhaj [l938D . This relation 
is widely applied in astrophysics and assumed to be robust, yet 
the empirical evidence in support of it is relatively weak (e.g. 
IProvencal et al.1 1 997l [2002 1 . This is because, despite it being com- 
paratively straightforward to constrain the radius of a white dwarf if 
its distance is known, it is difficult to reliably or independently de- 
termine its mass. For example, for a typical field white dwarf with a 
parallax based distance estimate, the mass can be determined from 
the radius and the surface gravity, following Equation [TJ where M 
is the mass, g the surface gravity, R the radius and G the gravita- 
tional constant. However, here the uncertainty in the mass estimate 
is at least twice the fractional error in the radius, resulting in signif- 
icant scatter amongst empirical data obtained in this way. 



M = gR 2 /G 



(1) 



Fortunately, for white dwarfs which are members of open 
clusters, both the distance and the radial velocity are constrained 
and thus it is possible to determine the gravitational redshift. The 
masses of these objects can then be estimated via Equation[2] 



M = -ufl/0.635, 



(2) 



where v is the gravitational redshift in kms -1 and M and R are 
mass and radius in solar units respectively. The uncertainty in the 
mass estimate is reduced since the fractional error in the radius de- 
termination has half the impact as in Equation[Tj 

In light of this, we have assumed that the candidate Praesepe 
white dwarfs which are consistent with the general trend of the 
IFMR, as defined by data points from other systems, are robust 
cluster members. In this case, the distance to each degenerate has 
been taken to be that to the centre of the Praesepe, but with an 
uncertainty associated with the intra-cluster spatial distribution of 
the white dwarfs. Here, we adopt D=184.5±6 pc as the distance 
to the cluster centre, which is the weighted mean of the Hipparco s 
based measurement, (m-M) =6.24±0.12 dMermilliod et alJll997D. 
the gro und based parallax measurement of dGatewood & de Jongd 
11994) . (m-M)o=6.42±0.33 and a recent photometric determination 
which takes into account a new spectros copic determin ation of the 
cluster metallicity, (m-M) =6.33±0.04 dAn et alj|2007h . All white 
dwarf candidate m embers of Praesepe identified in the survey of 
IDobbie et aHd2004l) . which extended beyond 2.5° of the cluster, are 
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Table 3. ID, name, coordinates, expected line-of-sight velocity assuming cluster membership (RV"c), weighted mean line core shift, gravitational redshift 
based on UVES T e g and log g, line-of-sight velocity (RVq) and estimated deviation from cluster velocity ((RVq-RVc)/(ARVo 2 + 0.8 2 ) ' 5 ), for the eight 
non-magnetic white dwarfs included in this work. 



ID 


Name 


RA 


Dec 


RV C 


H-a,H-/3 shift 


GR. 


RV ±ARV 


(RV -RV C )/(ARV 2 +0.8 2 ) 5 






J2000.0 






kms 






WD0833+194 




08:36:10.01 


+19:38:19.1 


35.0 


78.6±1.2 


40.1 ±4.4 


38.5±4.6 


+0.74 


WD0836+199 


LB1847.EG60 


08:39:47.20 


+19:46:12.1 


34.5 


86.1=1=1.7 


43.4±4.7 


42.7±5.0 


+1.61 


WD0837+185 


LB5959 


08:40:13.30 


+18:43:26.4 


34.8 


84.3±1.5 


49.2±5.2 


35.1±5.4 


+0.05 


WD0837+199 


LB393, EG61 


08:40:28.09 


+19:43:34.8 


34.5 


73.5±1.1 


41.5±4.5 


32.0±4.6 


-0.53 


WD0837+218 




08:40:31.47 


+21:40:43.1 


33.9 


81.9=1=1.1 


62.8±6.5 


19.1=1=6.6 


-2.24 


WD0840+190 




08:42:58.03 


+18:54:35.5 


34.4 


84.5±2.2 


54.7±5.7 


29.8=1=6.1 


-0.74 


WD0840+200 


LB 1876 


08:42:52.32 


+19:51:11.3 


34.2 


74.9=1=1.1 


40.1±4.4 


34.8±4.5 


+0.13 


WD0843+184 


LB 8648 


08:46:01.91 


+18:30:48.5 


34.2 


88.0=1=1.2 


51.5=1=5.4 


36.5±5.5 


+0.41 




0.4 0.6 0.8 1 

Mass (M Q ) 



Figure 6. The mass-radius relation for white dwarfs. The locations of the 
five individual Praesepe white dwarfs for which we were able to estimate 
radius and mass are shown (light-grey points) together with their weighted 
mean (large black circle). Three degenerate components of visual binary 
systems for which the most robust estimates of mass and radius are cur- 
rently available (black squares) and theoretical mass-radius relations for 
white dwarfs with CO cores with thick H layers at effective temperatures 
of T cff =5000 K (bottom), 15000 K (middle) and 30000 K (top) are also 
shown (solid black lines). 



found to lie within 2° of the cluster centre, corresponding to 6.5 pc 
at D=184.5 pc. If the distribution of the white dwarfs is spherically 
symmetric about the cluster centre then all are likely to lie within 
the distance range D=184.5±8.5 pc. 

The radius of each degenerate has been determined by scal- 
ing the distance (D) by the square root of the estimated flux ra- 
tio, (f/F) 0,5 , where f is the observed flux at the Earth's surface and 
F is the flux at the surface of the white dwarf. The surface flux 
in the Sloan Digita l Sky Survey (SDSS) filters (r, i and z; e.g. 
Fukug itaet al J 1996b for each star has been obtained f rom the pure- 
H model atmospheres of (Holberg & Bergeron 2006), where cubic 
splines have been used to interpolate between points in the grid. 
The uncertainty in the surface flux due to the error in the surface 
gravity determination for each object was calculated to be 0.3-0.4 



per cent, assuming a generous uncertainty of ±0. 1 in the log g mea- 
surement, and so was assumed negligible. However, in the temper- 
ature range of interest, the error in the effective temperature deter- 
mination for each object was found to introduce an uncertainty of 
~3.5 per cent in F (r, i and z). We note that this uncertainty would 
have been larger if the u or g bands were used. The observed SDSS 
fluxes (Table [4) are believed to have absolute uncertainties at the 2 
per cent level flAdelman-McCarthv et alj2008l) . Thus the estimated 
total combined uncertainty in our white dwarf radius determina- 
tions is ~5 per cent, where contributions from errors in the white 
dwarf distances, the predicted fluxes and the observed fluxes are 
4.6 per cent, 1.75 per cent and 1 per cent respectively and have 
been added quadratically. 

Subsequently, the mass for each star has been re-calculated 
using Equation [2] this time using a gravitational redshift deter- 
mined by subtracting the predicted line-of-sight velocity (RVc in 
Table |3), assuming each white dwarf shares the cluster motion, 
from the weighted mean measured shift of the non-LTE cores of 
the Ha and Hf3 lines (Table |4j- As open cluster escape velocities 
are typically low, app roximately twice the velocity dispersion (e.g. 
iFellhauer et alj2003h . any object that has a space velocity substan- 
tially different from the cluster mean will be rapidly lost from the 
system. We note in this context that as the white dwarf members of 
Praesepe have cooling ages of r coo i-~ 200-300 Myr, it is unlikely 
that putative point asymmetries in mass loss during the later stages 
of AGB evolution have displaced their space velocities greatly from 
the cluster mean. Moreover, th e velocity dispersion of ~1 M q clus- 
ter members as determined by Mermilliod & Mavorl |l999) is only 
0.8 kms -1 . Propagating these uncertainties, we estimate that the 
mass determination for each white dwarf using the gravitational 
redshift should be robust to ~6 per cent. 

The magnetic white dwarf, WD0836+201, the possible 
double-degenerate system, WD0837+185, WD0836+199 which 
has no clean SDSS photometry and WD0837+218 which sits well 
above the bulk of objects in the IFMR have been omitted from these 
calculations. Our estimates of the masses and radii of the remaining 
five white dwarfs are given in Table [4] 

These masses and radii are plotted in Figure [6] (light grey 
points with error bars) together with the most robust points cur- 
rently available in the literature (small points with error bars). 
The masses of Procyon B (T cff Ri7700 K), 40 Eri B (T off « 16700 
K) and Sirius B (T o fr~25000 K) have been determined precisely 
using dynamical methods since these white dwarfs are members 
of visual binaries. Moreover, the distances to these three systems 
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Table 4. The SDSS magnitudes, the model independent masses and radii for the white dwarf members of Praesepe with UVES spectroscopy. The magnetic 
white dwarf, WD0836+201, the possible double-degenerate system, WD0837+185, WD0836+199 which has no clean SDSS photometry and WD0837+218, 
which sits above the bulk of white dwarfs in the IFMR, have been omitted from these calculations. 



Name 


r 




i 




z 


M 


R 














M Q 


xlO" 2 R 


WD0833+194 


18.09±0.02 


If 


,.34±0.02 


18 


.65±0.04 


0.78±0.05 


1.133±0.057 


WD0837+199 


17.81±0.01 


If 


>.O8±0.02 


18 


.34±0.04 


0.71 ±0.04 


1.160±0.059 


WD0840+190 


18.44±0.01 


If 


>.70±0.02 


18 


.98±0.04 


0.75±0.05 


0.954±0.049 


WD0840+200 


18.01±0.01 


If 


,.21±0.02 


18 


.57±0.05 


0.75±0.05 


1.175±0.060 


WD0843+184 


18.36±0.01 


If 


,.61±0.02 


18 


.85±0.04 


0.84±0.05 


0.993±0.050 



are w ell constrained by Hipparcos parallax measurements (ESA 
1 19971) . Theoretical mass-radius relations for T cff =5000, 15000 and 
30000 K, based on the CO core, thick H-la yer evolutionary mod- 
els of Fontaine. Brassard & Berg eror]<r200lh are also shown (black 
lines). Four out of the five Praesepe stars are located within la of 
the T cff =15000 K track. WD0840+190 lies somewhat below this 
(~2er). It is possible that the radius of this object has been un- 
derestimated perhaps because it resides on the very back edge of 
the cluster. However, there is good agreement between the spectro- 
scopic and flux based radius estimates for all five white dwarfs. Al- 
ternatively, the radial velocity of this object may be overestimated 
by ~4-5 kms - , resulting in a marginally low gravitational red- 
shift determination. WD0840+190 may have very recently received 
a velocity kick through an interaction with a binary system near the 
centre of the cluster and be in the process of leaving Praesepe. A 
remaining possibility is that the low-level ripple in some of our 
UVES spectra could have adversely affected our determination of 
the shift of the line cores. Given that the slope of the IFMR is not 
expected to change significantly within the comparatively narrow 
mass window in which our objects reside, the impact of a number of 
uncertainties e.g. the velocity dispersion, the intra-cluster distribu- 
tion of the white dwarfs and the random errors in our line core shift 
measurements, can be reduced if we take the weighted mean of the 
masses and the radii of the five stars. In this case, we determine the 
M mcan = 0.760±0.021M Q andi? mcan = 0.01067±0. 00024 R Q 
(large filled circle in Figure^. 

Given that the mean effective temperature of these five objects 
is T e g~ 15400 K we find the location of this point in mass-radius 
space to be entirely consistent with the theoretical relation. Indeed, 
the proximity of all four robust data points to the model tracks, 
strongly supports the theoretical white dwarf mass-radius relation 
and confirms the veracity of the prediction that white dwarfs with 
larger masses have smaller radii. 



the mean cluster velocity due to evolutionary effects (e.g. asym- 
metric mass loss) since they have remained within the projected 
tidal radius of Praesepe despite having formed 200-300 Myr ago. 
While the gravitational redshift component to the velocity shift of 
the absorption line cores in the spectra of white dwarfs is apprecia- 
ble, given the robustness of evolutionary models as demonstrated 
above, if the effective temperature and surface gravity of a star is 
known, then the magnitude of this effect can be calculated straight- 
forwardly using Equation|2] 

Hence the radial velocity for each candidate Praesepe de- 
generate has been obtained by taking the difference between 
the measured shift of the line cores and the gravitational red- 
shift as determined from evolutionary models and the measured 
gravity and effective temperature. It can be seen from Table [3] 
that, as expected, the radial velocities of the bulk of the white 
dwarfs, including the mean velocity of our candidate double- 
degenerate system, lie within la of that expected if they are mem- 
bers of Praesepe. These latter values have been estimated as- 
suming a total cluster space velocity of Vtot=47.5 kins -1 , a 
conv ergent point of RA=06hl4m00s, Dec=-04d36mOOs J2000.0 
(e.g. IMermilliod et al.lll99fjl: IMermilliod & Mayor! 1 19991) and the 
method detailed in iReid Jl996h ." We find that the measured ra- 
dial velocity of only one star, WD0837+218, is notably different 
(14.8±6.6 kms -1 or 2.24a below), from the velocity predicted on 
the assumption of cluster membership. There is no evidence in the 
available data that this object is a radial velocity variable. As this 
white dwarf has the largest projected separation from the cluster 
centre of those identified by Dobbie et al. (2004), and appears to 
be an outlier in the IFMR of iDobbie et alj (2006a) we are led to 
conclude that it is more likely a field star interloper than a member 
of Praesepe. Assuming this to be correct, it should not be used for 
constraining the form of the IFMR. 



3.4 Radial velocities and the membership status of 
WD0837+218 

As has been discussed in the introduction, the white dwarf mem- 
bers of open clusters are extremely useful for constraining the 
form of the IFMR. However to assess the level of intrinsic scat- 
ter in this relation it is crucial that field white dwarf interlopers 
are eliminated from samples. An excellent way to adjudge clus- 
ter membership is via the measurement of radial velocities. Open 
star clusters typically hav e smal l velocity dispersions. For exam- 
ple, [MermTlliod & Mavor (1999) determine a velocity dispersion 
of only 0.8 kms" 1 for a sample of single F5-K0 members of Prae- 
sepe. Moreover, as discussed previously, it is unlikely that these 
white dwarfs will have velocities which differ substantially from 



3.5 Minor revision to the Initial Mass-Final Mass Relation 

These new UVES data for eight white dwarf members of Prae- 
sepe are only of moderate S/N per pixel but are homogeneous 
in nature and of much higher resolution than previously pub- 
lished spectra. The line shift measurements suggest that previous 
studies of these stars may have systematically underestimated by 
a small amount (~0.1) their surface gravities. For example, for 
WD0843+184 we derive a radial velocity of 45.2 kms -1 based 
on the surface gravity (and effective temperature) of Dobb ie et al.l 
(2006a). Additionally, wc estimate a radial velocity of 42.9 kms -1 
for WD0837+199 adopting the surface gravity (and e ffective tem- 
peratu re) erroneously assigned to WD0836+201 in Clav er et"aT1 
(2001). On this basis we argue that our measurements of effective 
temperature and surface gravity are, except in the case of the mag- 
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netic white dwarf WD0836+201, to be preferred to those obtained 
previously. 

As in our earlier work (e.g. iDobbie et alj|200 6a). we have 
utilised a grid of evolutionary models based on a mixed CO core 
composition and a thick H surface layer (e.g. lFontaine et al"]|200ll) 
to estimate the mass and cooling time of each white dwarf from 
our measurements of effective temperature and surface gravity (see 
Table QJ. Cubic splines have been used to interpolate between the 
points within this grid. The lifetime of the progenitor star of each 
white dwarf has then been calculated by subtracting the cooling 
time from the age of the cluster (625±50 Myr). Subsequently, to 
constrain the mass of the progenitor star, a s in our earlier wor k 
we have used the stellar evolution models o f iGirardi et al] £2000), 
here for sol ar metallicity (Zp)=0.019, which is close to the value 
favoured by lAnders & Grev esse 1989). Again cubic splines were 
used to interpolate between the points in the grid. We find that these 
model grids indicate that the difference between the lifetimes of 3-4 
Mq stars with Z=0.019 and Z=0.030 compositions is only ~1 per 
cent. As the metallicity of Praesepe is midway between these val- 
ues, the use of calculations for solar composition to estimate initial 
masses here appears adequate. Our estimates of the masses, cool- 
ing times and progenitor masses for each of the white dwarfs are 
shown in TableQ] 

The revised locations of the Praesepe white dwarfs are 
plotted in Figure [7] where the cooling times, masses and ini- 
tial masses for WD0836+197 (LB5893), WD0840+205 and 
WD0836+201 are based on effective te mperature and sur- 
face gravity estima tes from the literature dClaver et al. I l200ll ; 
IDobbie et all l2006ah . WD0837+218 has been excluded since 
based on the arguments presented above it appears more likely 
to be a field star. We also include in this plot the hydrogen 
rich white dwarf members of a number of other systems where, 
for a given initial mass range, parameters such as cluster 
age are re latively well const rained - namely NGC6819 and 



As has been shown by a number of authors (e.g.lFerrario et al 



NGC7789 1 Kali rai et al] |2008|). NG C6791 dKalirai et al] 120071). 
NGC 6633 (Willi ams & Boltel 120071). the Hyades teergeron etatl 



19951) . NGC2168 (M35: I Williams . Bolte & Koesten l2004 
2008). NGC2516 dKoester & Reimersl 1 19961) and the Pleiades 
( iBergeron et alj [l995; Dobb ie et alj|2006al lbl)~and the Sirius binary 
dLiebert et al]|2005l) . For the degenerate members of the Hyades, 
NGC2168, NGC2516, the Pleiades and the Sirius binary we have 
made the same assumptions regarding white dwarf effective tem- 
peratures and surface gravities a nd cluster ages and metallicities as 
we did in lDobbie etal] d2006al lbh. For NGC6819, NGC7789 and 
NGC6633 we have assumed solar metallicity and have adopted the 
cluster ages and white dwarf effective temperatures and surface 
gravities given in the referenced literature. We have excluded the 
three lower mass white dwarf candidate memb ers of NGC1039 
dRubin etal]|2008l) and NGC7063 LAWDS 1 dWilliams & Boltel 
120071) since pr oper motion measur ements indicate these are proba- 
bly field stars dDobbie et al1l2009t) . For the super-sol ar metallicity 
NGC6 791 we have used the [Fe/H]=+0.18 models of IGirardi et al] 
(2000) since this is the most metal rich composition available in 
their grid. We have included only the white dwarfs with a mass 
indicative of a CO core. The lower mass helium white dwarfs 
in th is cluster could be the products of close binary evolution 
(e.g. iBedin et al] [20081 : | van Loon et al. I 120081) or the products 
of a metallicity enh anced mass l oss me chanism dKalirai et al] 
120071) . Note that as in lDobbie etal] d2006al lbl). the NGC2168 white 
dwarf candida t e mem ber LAWDS 22 was not included here as 
Willi ams et ail d2004l) were unable to obtain a satisfactory model 
fit to the observed Balmer lines. 



2005 ; IDobbie et al] l2006al ; Iwilliams & Boltel [2007 ; Kalir ai et al 
2008), the IFMR as defined by the vast majority of current semi- 
empirical data (Figure[7) can be reasonably approximated by a sim- 
ple linear function (dashed line). Interestingly, WD0837+185, the 
radial velocity variable, does not appear to be an outlier. If the sys- 
tem experienced a common envelope phase then it may well have 
occurred during the later stages of the AGB. In this case the evolu- 
tion would approximate that of a single star. Indeed, there is only 
one star, WD0836+197, which deviates substantially from the gen- 
eral trend. Intriguingly, there appears to be nothin g else untoward 
about this star. For example, Dob bie et al] |2006a) confirmed it to 
be a pr oper motion member, while the spectroscopy of lClaver et al] 
reveals no evidence of Zeeman shifts in the Balmer line 
cores. The line-of-sight vel ocity as der ived from the line core ve- 
locity shift measurement of Reid] d 19961) and the spec troscopic ef- 
fective temperature and surface gravity estimates of Idaver et al] 
d200ll) is extremely close to the value expected on the basis of clus- 
ter membership. This argues against, b ut does not rule o ut, it being 
a radial velocity variable. Furthermore. iKarl et al] 120051) found no 
evidence for rapid rotation, the shape of the Ha line core setting a 
limit of v sin i <22 kms -1 . 

The magnetic white dwarf, WD0836+201, sits marginally 
above the bulk of white dwarfs in Figure UJ However, T B s and 
log g from line fitting are likely to be less robust here than for 
non-magnetic degenerates. Indeed, the SDSS photometry (g — i, 
u — r) is more consistent with a slightly lower effective temper- 
ature (T c ff~ 15000 K). We note that despite a gravitational red- 
shift based mass estimate ( Heber, Napiwotzki & Reid 1 19971) ap- 
pearing to_be_randstentwith the spectroscopic mass determina- 
tion of Idaver et al. I feooib . this result is void due to the confu- 
sion between WD0836+201 and WD0837+199. In any case, it is 
not obvious that a strongly magnetic white dwarf should be ex- 
pected to follow an IFMR delineated by non-magnetic stars. Mag- 
netis m may hinder mass loss during post-main-sequence evolu- 
tion dWickramasinghe & Ferrario 2000). Alternatively, all strongly 
magnetic white dwarfs ( B<;1 MG) might be produced through 
close binary interaction. iTout et ah |2008) propose that isolated 
high field magnetic white dwarfs form via the merging of two stel- 
lar cores within a common envelope environment. This unusual 
evolutionary history would presumably result in a higher than ex- 
pected white dwarf mass for the estimated progenitor mass. 

Thus the location of WD0836+197 in initial mass-final mass 
space may represent the best evidence for strong differential mass 
loss. However the lack of other comparatively deviant points in Fig- 
ure UJ and the sharply peaked form of the white dwarf mass distri- 
bution leads us to still prefer an evolutionary scenario in which it 
has formed from a blue straggler star. Praesepe is known to contain 
a number of blue straggler s tars (e.g. 40 Cancri and Epsilon Cancri; 
lAhumada & Lapassetll2007l) . 

A revised linear least squares fit to the 41 white dwarfs in 
this plot, which excludes the magnetic white dwarf WD0836+201 
and the outlier WD0836+197, between progenitor masses of 
1.15 M and 7 Mq gives parameters m=0.1197±0.0074 and 
c=0.3569±0.0220. These parameters are not greatly different from 
our original estimate s of m=0.133±0.015 and c=0.289±0.051 
dDobbie et ail l2006ah or those from the recent work of other 



groups e.g 
c=0.394±0.025. 



Kalirai etal] d2008l) determine m=0.109±0.007 and 
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Figure 7. A plot showing the revised locations of the Praesepe white dwarfs in initial mass-final mass space (filled circles). The locations of the degenerate 
members of the NGC6633 (square star), NGC7789 (cross), NGC6819 (asterisk), Hyades (open triangles), Sirius binary system (open circle), the Pleiades 
(open stars), NGC2168 (open diamonds), NGC1039 (open squares), NGC6791 (filled triangles) and NGC2516 (open '+' signs) are also shown. The relation of 
Weidemann 12000) (dotted line) and a revised linear fit to 41 white dwarfs (dashed line) are overplotted. The radial velocity variable, WD0837+185 is labelled. 



4 SUMMARY 

We have presented an analysis of high resolution VLT and UVES 
spectroscopy of nine white dwarf candidate members of Praesepe. 
For the eight non-magnetic stars we have measured the velocity 
shift of the Ha and H/3 line cores and obtained revised estimates of 
their effective temperatures and surface gravities. We have shown 
that contrary to the conclusions of a number of previous studies, 
WD0836+201 is the magnetic white dwarf. We have also provided 
evidence that WD0837+185 is a radial velocity variable and pos- 
sibly a double-degenerate system. Additionally, we have used the 
five most appropriate Praesepe stars to obtain a new robust con- 
straint on the white dwarf mass-radius relation. This new data point 
and the three most robust existing data points have been shown to 
lie in close proximity to the theoretical tracks. This result confirms 
the veracity of modern theoretical white dwarfs mass-radius rela- 
tions. In light of this we have determined the radial velocities of 
all eight non-magnetic degenerates. We have found that the line- 
of-sight velocity of WD0837+218, when examined in conjunction 
with its projected position with respect to the cluster centre and 
its location in initial mass-final mass space, argues that it is more 
likely to be a field star than a cluster member. It should thus be 
discounted from investigations concerned with the form of the ini- 
tial mass-final mass relation. We have employed our new param- 
eters for the Praesepe white dwarf members, in conjunction with 
comparatively robust data from other open clusters to re-assess the 
form of the initial mass-final mass relation. We find that there is 
only one substantial outlier, WD0836+197, from a near monotonic 
IFMR, which can still be adequately represented by simple linear 
function. The location of this star in initial mass-final mass space 
may be evidence of strong differential mass loss but given the lack 
of other outlying stars and the sharply peaked nature of the white 



dwarf mass distribution, we still favour an evolutionary scenario in 
which it has formed from a blue straggler star. 
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